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Abstract 



We present a new variational approach to the study of phase tran- 
sitions in frustrated 2D XY models. In the spirit of Villain's approach 
for the ferromagnetic case we divide thermal excitations into a low 
temperature long wavelength part (LW) and a high temperature short 
wavelength part (SW). In the present work we mainly deal with LW 
excitations and we explicitly consider the cases of the fully frustrated 
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triangular (FFTXY) and square ( FFSQXY) XY models. The novel 
aspect of our method is that it preserves the coupling between phase 
(spin angles) and chiral degrees of freedom. LW fluctuations consist 
of coupled phase and chiral excitations. As a result, we find that for 
frustrated systems the effective interactions between phase variables 
is long range and oscillatory in contrast to the unfrustrated problem. 
Using Monte Carlo (MC) simulations we show that our analytical cal- 
culations produce accurate results at all temperature T\ this is seen at 
low T in the spin wave stiffness constant and in the staggered chirality; 
this is also the case near T^. transitions are driven by the SW part 
associated with domain walls and vortices, but the coupling between 
phase and chiral variables is still relevant in the critical region. In 
that regime our analytical results yield the correct T dependence for 
bare couplings (given by the LW fluctuations) such as the Coulomb 
gas temperature Tqg of the frustrated XY models . In particular we 
find that Tcg tracks chiral rather than phase fluctuations. Our results 
provides support for a single phase transition scenario in the FFTXY 
and FFSQXY models. 

PACS Numbers 75.10-b, 75.10.Hk 
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INTRODUCTION 



Frustrated magnetic systems have been extensively studied, in part because 
they constitute non-disordered versions of spin-glassesBH. They display rich low- 
temperature phases and remarkable phase transitions since frustration modifies the 
naive symmetry of the Hamiltonian. For spatial dimensions D > 3, Kawamura 
found by renormalization group (RG) techniques that frustrated 0{n) spin mod- 
els belong to a new, "chiral" universality clasdl. For D = 2 and XY spins, which 
is our present concern, phase transitions are dominated by defects. Frustration 
results in additional chiral variables, which generate a discrete symmetry. In the 
fully frustrated case the Hamiltonian for a square lattice is believed to possess an 
0(2) X Z2 symmetry (see for instance Reli); for the triangular lattice (FFTXY) one 
has the extra Csv symmetry associated with the permutation of the three sublat- 
tices, thus adding the possibility of a Potts transition!. The transition associated 
with the 0(2) part (phases i.e. angular variables) would be Kosterlitz-Thouless 
(K-T)-like at a temperature T^-tJ^ and the discrete Z2 part (chiral variables) would 
be broken below a temperature T^s- There is an ongoing controversy concerning 
the order in which these transitions should take place. RG calculations suggest 
that Tq = Tkt = Tds but two transitions are not ruled outBi"0 : in the single 
phase transition scenario, measurements of critical exponents for the chiral and of 
the central charge using MC and MC transfer matrix techniques^ reveal non Ising 
behavior, providing support for Kawamura's claim of a new universality class even 
in 2D. This is also suggested by studies based on the selective breaking of certain 
symmetries£3ll3. Some MC studies performed on the FFTxyS'ilil and on the fully 
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frustrated square (FFSQXY) mo yield a single phase transition; yet other MC 
simulations for the FFSQXY model and for the 1/2 integer Coulomb gas give two 
phase transitions very close in temperatureEH^l. 

In view of these unsettled issues the present paper has two objectives : 

First we would like to give a quantitative description of the relevant excitations 
in FF systems. In doing so, we wish to assess the importance of the coupling between 
phase and chiral degrees of freedom at low temperature (T) and in the critical region. 
This would allow us to identify the nature of the critical fluctuations and to decide 
whether one should expect two phase transitions or just one. 

Second we would like to test if the thermodynamic properties of the FFTXY 
and of the FFSQXY models are similar or not, and in particular if the nature of the 
phase transitions is different or the same for the two systems. 

In order to get some insight into these issues we present a new analytical ap- 
proach to study 2D XY frustrated systems. It is inspired by Villain's analysis for 
ferromagnetic (F) systemsii where: (a) the partition and correlation functions are 
products of a Long Wavelength (LW) - spin waves - and of a Short Wavelength 
(SW) - vortex - contribution@; and (b) the long wavelength part is mapped quan- 
titatively (by perturbation theory or variational scheme) onto the low temperature 
contribution of the original cosine Hamiltonianii. Steps (a) and (b) allow one to 
compute accurately all relevant thermodynamic quantities from T = up to T^. 
Our results are best summarized on the figures. In section I we begin with a brief 
discussion of steps (a) and (b) for the unfrustrated case and we show that a sim- 
ple variational approach - the self-consistent harmonic approximation (SCHA)i - 
yields quantitative agreement with MC at all T so long as one considers thermody- 
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namic variables sensitive to LW excitations (Fig (1)). Extending the method to FF 
systems produces incorrect results even at very low T (Fig (2)). Failure is due to 
the fact that a naive application of the variational approach eliminates chiral fluctu- 
ations. In section II we set up a new variational method ( we call it NSCHA for new 
SCHA ) which explicitly preserves the coupling between phase and chiral degrees of 
freedom. As a result LW excitations consist of coupled spin waves (phase) and polar 
(chiral) fluctuations. This feature causes the effective interactions between phase 
variables to be long range and to oscillate in sign (Figs (3)- (4)). One may contrast 
this behavior with the unfrustrated case where phase couplings remain short range 
and positive in sign. Focusing on LW fluctuations, in section III we compare our 
results to MC simulations performed on the FFTXY and FFSQXY systems. For 
the FFTXY model we have done MC calculations for sizes up to 60 x 60, using 
10^ - 10^ MCS/spin and fluctuating boundary conditionsElS. For the FFSQXY 
lattice we have used available data from the literature : this is justified since we 
only use in a quantitative fashion data pertaining to LW excitations (i.e. results 
which are common to all studies ). For all T MC simulations and analytical results 
agree closely: this is apparent in the LW contributions to the stiffness constant and 
to the chiral order parameter (Figs (5)- (8)). Above a characteristic temperature T* 
defects become important and they ultimately drive the transition: we introduce a 
variable r equ ( p6D which allows to track chiral domain walls. Fig (7) shows that 
they become relevant above T*. Furthermore Figs (5) and (6) show that for T > T* 
domains affect both the chiral order parameter and the spinwave stiffness constant. 
The relevance of the coupling between phase and chiral variables near the transitions 
is also visible on Fig (8) where we can see the agreement between analytical and MC 
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predictions for the bare couplings - here the Coulomb gas temperature Tcg (equ 



37|) -. Noteworthy is the fact that Tcg is connected to chiral variables in FF systems 
whereas it is a bare coupling constant for phase fluctuations in the unfrustrated 
situation (Fig (1)); this point is important in view of the fact that MC studies of FF 
systems assume that Tcg is a bare coupling for the phase variables (see discussion 
in section III). Our study thus suggests that the vanishing of the spinwave stiffness 
and of the chiral order parameter occur at the same T, for the FFTXY and for the 
FFSQXY models. 



I. THE SELF CONSISTENT HARMONIC APPROXIMATION 
A. The ferromagnetic case 
Using standard notation the Hamiltonian reads 

^ = - E JijCosiOi - Oj) (1) 

<i,j> 

In the ferromagnetic case [Jij = J > for nearest neighbor pairs) Villain replaces 
the cosine potential by a parabolic form; ^PJ^^^i^i-^j) jg approximated by : 

Const. J2 e-'^-^^^^'-^^-^-"^^)' (2) 

where the integers Uij express the periodicity of the original interaction. The 
main features of the "Villain form " are that it includes both LW excitations (spin 
waves connected to the phases 9) and SW excitations (vortices connected to the lat- 
tice curl of the n) and that the partition function is the product of the LW and SW 
parts. Below Tkt the vortex part is essentially irrelevant ( it simply introduces a 
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dielectric constant ey ~ 1 + e~'^/"^) and LW properties are described by an harmonic 
spin wave hamiltonian with a spinwave stiffness constant 7 = Jy/ey. The impor- 
tance of Villain's form stems from the fact that by applying the Migdal-Kadanoff 
scheme to the original cosine interaction, one iterates towards an effective harmonic 
spin wave theory below T/^t^- The shortcoming of (2) is that the coupling con- 
stant of the LW part is temperature independent whereas the original cosine form 
introduces interactions between spin waves. To deal with this issue one may use 
the self consistent harmonic approximation (SCHA)§ : for the LW part one uses a 
variational hamiltonian 

Anharmonicites of the cosine potential translate into a temperature dependence 

of Jij. 

The variational free energy is given by F-^ar = Fo+ < — TYq >Ha (-^o is the free 
energy for hamiltonian Tio ) and reads 

Fvar = -IT. J^mJ ' E J^J(^~~''" + ^ log DctiJ/T) (4) 

<hj> <i,j> 

J is the matrix with diagonal elements ^^Jik and off diagonal elements —Jij- 
The quantities yij =< {6i ~ 6jY >Uo are themselves functions of the variational 
parameters Jij. Therefore we may use Uij as alternative variational parameters. 
The variational equations read : 

Jij = Jije'^y^' (5) 

with 
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^ [f ^ - cosq.{ri - fj)) 

Here J{q) is the Fourier transform of Jij. In the ferromagnetic case, the effective 
interactions Jij only couple nearest neighbors. Their magnitude J{T) is given by 



J(T) = Je (7) 

{z is the number of nearest neighbors of the lattice). One may then compute the 
spinwave stiffness matrix; its elements are given by the second derivatives of the free 
energy with respect to uniform twists of the phase0. For isotropic lattices the matrix 
is diagonal and all the elements are equal. In SCHA the constant is 'Jscha = J{T)- 
We may then simply replace Jy in equ (0) above by ■yscHA- Figs (la) and (lb) 
show the temperature dependence of the SCHA stiffness for the square (SQ) and 
triangular (TR) lattices along with the MC result (denoted by 7(T)). At low T 
vortices contribute with a probability of order e"^"^/^ so that SCHA and MC agree 
quite well. Near Tkt ~ 0.892J (SQ latticdlHS) or 1.446J (TR latticdi) vortices 
cause a drop in 7(T). Since SCHA only describes LW fluctuations it fails to produce 
a fall-off. 



B. the Fully Frustrated case 

The previous analysis is easily extended to the situation where spins are non 
collinear in equilibriumll. We rewrite the 6i of equ (|lD as: 

0^ = + ^^ (8) 



where 6i =< 6i >'Ho ^-^d the variational hamiltonian is 

^0 = ^ E (9) 

^ <i,i> 

In addition to the parameters yij =< {(pi — (pjY >Ho has the extra variables 
Oi' . The variational equations now read 

Jij = JijCosiOi^ - %°).e-5?/». (10) 
E Jijtanie^" - = (11) 

<i,j> 

and Hij is given by equ.(|^). Denoting by = [xi, Ui) the vector connecting the origin 
of the lattice to site i and by Uij the vector connecting nearest neighbor sites i and 
j, the solution OP to these equations is independent of T and given by: 

- e,' = e\n + u,,) - 9%^) = a,, + {-If^+y^asQ {mod 27r) (12) 

for the square lattice, along with the symmetry property 

e\n + u,,) - e\n) = +{e\n - u,,) - e\n)) (13) 

and by: 

^/ - = e\fi + Ui,) - e\n) = Q.Uij = Aj ± aTR {mod 2tt) (14) 

for the triangular lattice, along with the symmetry property 

^°(r-; + u,,) - 9\ri) = -{e\r, - Ui,) - e\n)) (15) 

Here Aij = —Aji = vr if Jij < and Aij = if Jij > 0. a is a lattice dependent, 
temperature independent quantity. For the FFTXY lattice {Jij = —J < ) one 
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has Q oc (f , so that utr = f . For the FFSQXY lattice - the so-called Villain 
odd model - (Jjj = — J < every other row along say the horizontal direction and 
Jij = +J > otherwise) one finds asQ = f • In addition Jij is a nearest neighbor 
interaction of magnitude 



J(T) = Jcos{a).e (16) 

Fig (2) shows the SCHA stiffnesses for the FFTXY lattice together with the MC 
result. Even at low T, J{T) and 7(r) differ significantly. The same effect is observed 
for the FFSQXY lattice. The reason for this discrepancy is clear : inserting equ (||) 
into equ (j^) gives 

= — ^ Jij{cos{Oi' — 9j^)cos{Lpi — Lfj) — sin{6i' — 6f)sin{Lpi — Lpj)) (17) 

<i,j> 

Within SCHA the sin() term of equ ([17D averages to zero. But this term precisely 



discriminates the +a solution from the —a solution in equs (|T2|,|T^) and these two 
solutions correspond to the two chiral groundstates of the FF system. As a result 
SCHA washes out chiral fluctuations and maps the hamiltonian onto an effective 
ferromagnetic phase problem since the Jij are simply renormalized to JijCos{6i^ — 

ef). 

Our analytical method was thus required to preserve the coupling between the 
two chiral states and to allow fluctuations of the chiralities. Next section shows 
that our approach then yields accurate results for the FF case. Furthermore it also 
improves on the standard SCHA in the unfrustrated case. 
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II. NSCHA 



Using equ (|T^ the partition function reads : 

Z = Tr^M^i - ^j)e''^<''^> J^j-^^^'-^^'y^^^^-^^) (18) 
where I2 = e~^^<^-^> )sin{ip^-Lpj) ^ rewrite I2 as the sum of a term even 

in ip plus a term odd in ip : 

h{x) = \[h{x) + h{-x)] + \[h{x) - h{-x)] 

Now the trace over p in equ ( pTSD is constrained by equ (H) Oi = Oi' + ipi but for 
LW excitations we expect pi to fluctuate about 0, so that we can safely extend the 
domain of variation of pi to the interval [— vr, +7r]. As a result the odd term of I2 
drops out and the partition function reads 

Z = Tr^o.e^^^^^'^' where 

{19) 

-TLog[cosh{J2<i,j> /3JijSin{9i° - 9j°)sin{(pi - ipj))] 
The second term on the r.h.s of equ (|l^) is the new relevant term. It has the 
following properties : 

• it is a thermal contribution since it vanishes at T = 

• it is zero for coUinear systems (i.e. either unfrustrated or frustrated but with 
non chiral configurations) 

• it couples chirality (6*° - more precisely Jijsin{9i^ — 9j^) for each link (ij) equ 
( pOD below - ) and phase {(pi) variables, allowing chiral fluctuations. Similarly it 
preserves the symmetry between the two chiral groundstates. 

We now compute the variational free energy associated with Tie// using the trial 
hamiltonian Hq (equ (|)). To do so we have to expand the log{cosh) term of Tig// 
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in power series of its argument. It is justified since the series amounts to a multipole 
expansion as is seen below. Besides, the leading term is the first term, especially at 
low T. 'Heff then becomes 

Ti-NSCHA = - T,<i,j> JijCOs{9i^ - 9j°)cOs{(fi - Lpj) 

The variational equations for - what we call - the NSCHA (new SCHA) ensemble 



are: 



cosh{yij + yki - yu - yjk) e"^(2'"'=+^^') 

Again yij =< (yj^ - (^jf and 

T ff J {1 — cosq.{ri — fj)) 



(21) 



:22) 



y = __ // d2q ' ---■^■l (23) 
{27t^)JJbz J(0)-J(g) ^ ' 

For the FFTXY and FFTSQXY lattices it is easy to check that Of is a temper- 
ature independent quantity and that its value is still given by equs (|T2|,|H) 



Futhermore equ ( PTj) shows that Jij is no longer a short range interaction. In 
fact we find that for all T, 

4- ~ 1^ ^^1. (24) 
12 



for large distances r = \ri — rj\ (Fig. 3). This comes about because i/ij ~ log{r) 
at large distances so that Uik + Vji — yu — Vjk ~ in equ (pTf ) : this contribution 
is quadrupolar like. Similarly, expanding the log(cosh) term to next order would 
produce a higher order multipolar contribution (see also Appendix A). 

In addition the sign of varies with the relative orientation of i and j and, in 
the case of the FFSQXY lattice, with the distance between i and j Fig (4a) and 
(4b). 

These features are to be contrasted with the results of SCHA yielding a positive 
nearest neighbor Jij. The coupling between phase and chiral degrees of freedom has 
produced an oscillating, "long range" interaction between the phase variables. Be- 
cause of these properties it is clear that renormalization group analyses (e.g Migdal- 
Kadanoff) are not straighforward for FF systems. 

Within NSCHA we can compute the phase stiffness constant. Owing to the 
isotropy of the lattices we have : 



(J(0) - J {q.u,)) 



r(T)= hm /'^^"^"^^"-"-^^ (25) 



where Ux is the unit vector along the horizontal direction of the lattice. Besides, 
within this new variational ensemble we also get a stiffness associated with the 
canting of the spins; considering a small variation of the nearest neighbor angle 
difference Oi' — Of' from its equilibrium value in the form Aux-Uij we get: 

(26) 

E<^,j> J:<k,l> J^JJkl (M,,-.M.)(uH.n,)e-2fe».+2^'='+?'«'=+?'.'-?^«-?^.'=) 

cos(^,° - %°)cos(4° - 61°) + siniOi'^ - ej^)sin{ek^ - 61°) 
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It is easy to show that '^nscha{T) is nothing but the average of the exact spin- 
wave stiffness 7(T)i in the ensemble Tio, i-e 

1nscha{T) = ^ < E<jj> Jij cos^Oi - 6j) {uij.Uccf 

-j^ 'E<i,j> E<fc,/> JijJki {uij.Ux){uki.Ux) sin{9i - 9j)sin{6k - 9i) >no 

(27) 

A plot of r(T) and 'jNSCHAiT) versus T is shown for the TR lattice (Fig (5a)) 
and for the SQ lattice (Fig (5b)). We note that r(T) and 'jNSCHAiT) coincide at 
low T. This is explicitly demonstrated in Appendix B. In particular we find that 

-fNSCHAiT) = 7o(l - ^) (28) 

• For the triangular lattice 70 = V3/2J and Tco = — J ~ 1.975J 

27r 

• For the square lattice 70 = \/2/2J and Tco = ~ 2. 075 J 
For both cases T^o ~ 2 J (ReS). 
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Another quantity of interest is the staggered chirality 

^ _ 1 < Eik,i)&pcrki >no ^29) 
A^p p E{k,i)&p (^kiiT = 0) 

P denotes plaquettes of the same sublattices i.e plaquettes in the same chiral state at 
T = 0. The summation J2{k,i)(^p is performed over the links of plaquette P oriented 
clockwise, aki is defined as : 

(^ki = Jkisin{9k - Oi) (30) 
(see below for a discussion on the definition of (Tm ). 
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Using equ ®, we have 



(^ki = Jki{cos{9k - 9i )sin{ipk - ipi) + sm(6'fc - Oi )cos{ipk - (pi)) 

Within NSCHA the sin{ifk — ^i) term drops out and a^iiT = 0) = Jki{cos{9k' — di^) 

so that 



where k and I are nearest neighbors {y^i has the same value for all the nearest 
neighbor sites / of any given site k). aNscHA versus T is plotted on Fig (6) for the 
TR and SQ lattices. 

To summarize the results of this section we see that LW thermal excitations in 
fully frustrated lattices are characterized by a strong coupling between chiral and 
phase degrees of freedom. The effective interaction between phase variables is long 
range and oscillatory - in contradistinction with the unfrustrated case -. 

Let us now compare our results to those coming from Monte Carlo simulations. 



In order to test the predictions of NSCHA we used the results of Monte Carlo 
simulations. For the FFSQXY lattice we took data from the literature insofar as we 
did not seek to extract information about critical fluctuations. In that case there 
is agreement among the various studies. For the FFTXY and for the ferromagnetic 



O' NSC HA — 




(31) 



III. MONTE CARLO VERSUS NSCHA 



triangular lattice recent data is rather scarce 



(Refs0'illi) so that we performed 



our 
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own simulations. We considered typical lattice sizes of 48 x 48 and ran 10^ — 10^ 
MCS/spin. In order to minimize boundary effects we used fluctuating boundary 
conditions (Refs0@). We monitored the following quantities : 

• the spinwave stiffness constant ■yiT) 

• the staggered chirality a equ (|29|) : 

several definitions of cr have been used in the literature. One of them is the definition 
we use here (see also Olssonll), others ar 

a^ = ^<J2Szgn{ M> (32) 

where aki is defined in equ (^) and Afp the number of plaquettes of each sublattice, 
or 

^' = ;^<E E <i> (33) 



with! 



P P {k,l)(^P 



a',i = ^{e, - 01 - A,i) (34) 

or with 

^H = ^iOk-ei) (35) 

In equs (0,^) the angular determination of the terms in parenthesis is taken in 
the interval ] — vr, +7r]. All these definitions lead to the same T dependence for a in 
the critical region. For the square lattice this is reported byEl for instance and for 
the triangular lattice this is seen on Fig (7) using the definitions equs (|30|) and (|35D . 

• the chirality amplitude r : 
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^_ 1 >p < Abs{E{k,i)ep(^ki) > .gg. 
Xp^ AbsiEik,i)epMT = 0)) ^ ' 

using again the previous definitions for (Jki ■ So long as cliiralities are ordered on eacfi 
sublattice r and a coincide. Wlien domains of the "wrong" chiral state form on a 
given sublattice the two quantities differ. Thus r allows us to track the formation of 
domains and domain walls (where the chirality of a plaquette J2{k,i)eP'^ki = 0). For 
the FFTXY lattice for instance Fig (7) shows that at Tc we have ~ 30% of positive 
chiralities, ~ 30% of negative chiralities and ~ 40% of a-chiral plaquettes on each 
sublattice. 

• the Coulomb gas temperature Tcg '■ 

This quantity monitors the bare (unrenormalized) coupling constant and allows 
to define the critical point for the XY model (Refs0S). Within MC it is given by0 

Tcg = TTT, Jo = J < cos{9i - > (37) 
for nearest-neighbors i and j. 



Figs (5a) and (5b) show T(T),ji\fscHA(T) and 7(T) versus T. The three curves 
yield the same variation at low T. Furthermore, as could be expected from our 
previous discussion, ■jNSCHAiT) tracks ■y(T) for T <T* (T* ~ 0.32J for the square 
lattice and T* ~ 0.35 J for the triangular lattice); for > j"* the two curves 
move apart. Since NSCHA describes LW excitations but neglects SW excitations 
responsible for the transitions, this had to be expected. 

Similarly Figs (6) shows a comparison between NSCHA and MC for o"(T); again 
the agreement is quite good for T < T*. Moreover we also see from Fig (7) that 
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T* marks the temperature above which domain walls become important, since cr(T) 
and t(T) start to differ for T ~ T*. 

At this stage we might worry that the discrepancy between MC and NSCHA 
predictions for T > T* not only marks the point when defects become important but 
also signals the breakdown of the variational approach. In fact NSCHA still yields 
accurate results for quantities sensitive to LW fluctuations including in the critical 
region. We see this by comparing the MC and the NSCHA Jq entering the definition 
of the Coulomb gas temperature. Tqg represents the bare (unrenormalized) coupling 
constant when LW fluctuations are taken into account. 

For instance in the case of the square lattice, Olsson finds that Tcg ~ 0.1280 
at the KT transition [T^t ~ 0.446 J); this value is to be compared with the MC 
results by GrestS {Tcg ~ 0.126) and by Ledsl (T^g ~ 0.1297) on the half integer 
Coulomb gas representation of the FFSQXY. 

Fig (8) shows Jq{T) for the square and triangular lattices determined both in 
MC and in NSCHA. We notice that: 

a) both determinations agree extremely well in the critical regime 

b) Jq tracks the chiral variable couplings rather than the phase variable couplings : 
Indeed, if we use for Jq the definition given in equ ( pT]) we find that in the NSCHA 
ensemble 

■JoNSCHA = cos{9^ - 6^)aj^sCHA (38) 

and Jq is therefore connected to the chiral variables. We have seen that the LW 
contribution to the chirality a - given by equ ( ^I]) - does not vanish at the transition 
(a becomes zero because of defects) so that Jqnscha is finite even in the critical 
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regime. For instance for the square lattice NSCHA gives Tqc ~ 0.125 using Tkt = 
0.446J. By contrast, for the unfrustrated case, equ (^) gives Jqscha = J{T) (see 
section I). So Jq is connected to phase variables then (see also below). 
In the frustrated case, if we replaced Jq by T(T) equ (|25|) we would find too high a 
value for Tcg (namely 0.139) compared to MC. 

Similarly, for the triangular lattice MC gives Tkt = 0.51J (Refsi'§) and Tcg ~ 
0.123 to be compared with the NSCHA prediction (using equ (|38D) Tcg ~ 0.122. 

This result has direct implications for MC studies : these introduce a second crit- 
ical temperature Tns where chiral order vanishes. Some authors find Tkt > Td^^ 
whereas others predict Tkt < TosSilii. At Tds critical exponents are found to 
be Ising-like by some authorslllS but non-Ising by otherffl. The magnitude of 
the jump of the spinwave stiffness constant or of the dielectric constant seen in MC 
appears universal for the FFTXY model (S) but non universal for the FFTSQXY 
modeSliUHll. Recently, Olsson has argued that a correct analysis of the transitions 
in the case of the FFTSQXY requires extra care due to their closeness in tempera- 
ture. As a result he finds a universal jump at Tkt and similarly Ising exponents at 
the chiral transitionil in contradistinction with previous author One should 

note that the claim of universality or non universality for the KT transition is based 
on a scaling a la Minnhagen for the magnitude of the jump : yet, according to 
Minnhagen's study this scaling should not hold (one might even expect a first order 
transition) given the value of the critical Coulomb gas temperature corresponding 
to Tkt^^- 

Our results show that, because of the coupling between phase and chiral degrees 
of freedom, Tcg pertains to chiral variables; because of this coupling one might thus 
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expect a single phase transition in these systems. 

For the unfrustrated case - e.g in the ferromagnetic hmit equ (^) gives 
Jqscha = IscHA = JiT). Its temperature dependance compares reasonably well 
with MC (Fig (1)). In fact MC and variational predictions agree extremely well if 
one compares Jo to 'Jnscha- NSCHA reduces to SCHA for the most part but even 
in the ferromagnetic case the stiffness 'Jnscha(T) equ (^) does not coincide with 
J(T). Chiral fluctuations exist even when 9^ = 0. Equ (|27|) shows that ■yNSCHAiT) 
represents the LW contribution to the stiffness constant 7, i.e the bare coupling 
for the phase variables. So it is natural to identify 2n-yi^scHA ^^^^ Tcg- If we use 
1nscha{T) in equ ( P7D we find analytically Tqg = 0.198 for the SQ lattice and 
Tcg = 0.191 for the TR lattice, to be compared with the MC values 0.1956 for the 
SQ latticed and 0.192 for the TR lattice. 

To summarize our results, we have constructed a variational ensemble (NSCHA) 
for fully frustrated XY systems in 2D. Testing its predictions with Monte Carlo 
simulations we see that our approach yields accurate results at all temperature - 
including in the critical regime - for quantities sensitive to long wavelength excita- 
tions. The key ingredient of the theory is the coupling between phase and chiral 
degrees of freedom and this coupling is always relevant. In particular it causes the 
interaction between phase variables to be polar-like (long range and oscillatory). 
As a result renormalization schemes assuming short range couplings might not be 
reliable. 

If a Coulomb gas temperature is introduced it appears to track chiral variables rather 
than phase variables. 

Monte Carlo simulations show that defects drive the transitions. In particular chiral 
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domains appear to affect tlie spinwave stiffness constant and ciiiralities in a similar 
fashion giving support for a single phase transition scenario. 

The above results pertain to both the FFTXY and the FFSQXY lattices suggesting 
universality for fully frustrated systems. 

For ferromagnetic systems NSCHA still improve on SCHA. The reason is because 
NSCHA incorporates fluctuations of the macroscopic phase (^o) about its equibrium 
(zero) value, in contradistinction with SCHA. In that sense NSCHA is a canonical 
ensemble as opposed to SCHA which is a microcanonical ensemble. In that limit 
the Coulomb gas temperature is associated with the bare coupling constant of the 
phase variables. 
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APPENDIX A: 



In this appendix we show that in the hmit of large distances R = \ fi — rj\, Jij 
given by equ (0) behaves as Jij ~ (equ (^). 

We start with the NSCHA variational equation for J^j : 



2T 



k,l 

1 



+ - J2 -^ikJji sin(6',° - el) sin(^J - O^) cosh(?/ij + - yu - yjk)e 

^ k,l 

where yij = y{rj — rj) is given by equ ( p3D and where the angles {0^} satisfy equs 



(0 0)- 

For the FFTXY the expression Jij sin(^^ — 9^) only depends upon Vj — . For the 
FFSQXY lattice however, there are four different types of sites (see equ (pT^)) so 
that Jjj sin(^° — 0^) exphcitely depends upon site i. Yet, for the FFSQXY lattice 
the quantity (— 1)^'+*'^ Jj^ sin(6'° — 6^) is independent of i . 



Therefore we introduce 



J,^.,^ = w^^+y'J,, sin(^° - 9^), (Al) 

For the triangular lattice we choose w = 1 and for the square lattice we set w = — 1. 
The quantity J defined in that way only depends upon fj — fi. J possesses the 
following symmetry properties : 

J(-e-) = -m (A2) 

for the TR lattice and 
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J{-e) = J{e), J{u,) = -J{uy) (A3) 

for the SQ lattice. 

In the following we introduce the notations R = fj — fi, e' = fi — fj , e = fk — Vi. 
For large R, only the third term in the r.h.s of equ ( pTD contributes, since Jij is a 
nearest neighbor interaction. Using the fact that Jij is independent of i, that we 
have the symmetry properties, equs ( [A2| , [A3D and that the value of yij is the same 



for all the nearest neighbors j of a given site i (by symmetry) we get 

Jn = -ie-^(^")^/;^^+^- Yl JeJp cosh \{y{R + e + e') + y{R) - y{R + e) - y{R + P)) 

where io is any nearest neighbor vector connecting two sites. 
We denote by Ay{R,e,e') the following quantity: 

Ay{R, e, e') = y{R + e + e') + y{R) - y{R + e) - y{R + e'). In the large R limit 
we expand y in powers of R; with : 

a, /J 

and 

-> 1 d^yi^R) 

^3y{R, ^' ^) = g X! (~^a^/3^7 ~ ^a^7^/3 ~ ^'(^^'■y^a + + ^a^-f^'iS + ^/3^7^a) np 

a,(3,y a /3 7 

we have that Ay{R, e, e*) = A2y{R, e, e') + Asy{R, e,e') + --- 



In Appendix B we show that Jo — Jg ~ Fg^ for small q so that - using equ 
we get y{R) log |^| for large R. As a result 

d^yjR) _1_ d^yjR) _1_ 

dRadRi3 ^ W dR^dRpdR^ ^ 'R? 
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(A4) 



Expanding cosh() we get: 



e',e 



that is 



(A5) 



Jr = -le-^(^°)^^^+^« Y: -Je-W + ll^^yiR, e, (A6) 
+\[A2y{R, e, P)] [A,y{R, e, e')] + l[Asy{R, e, P)]^} + 0{R-') 



• for the triangular lattice, using equ (|A2|) we have 

Ep^,JeJp[A2y{R,e,^')f = 0, 
Ep,,JeJp[A2y{R,e,e')][A3y{R,e,e')] = 0. 
so that 



(A7) 



• for the square lattice, using equ (^) we have 

Ep,,JeJp = 0, 

but Ep^,^JeJp[A2yiR,e,e')]^ ^ 0. 
so that 

We note the sign alternation due to the (—1)^^+-^!' term for the FFSQXY lattice 
(see Fig. 4) 
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APPENDIX B: 



In this appendix we compute ^NSCHAiT) (equ (p6D) and TiT) (equ (^)), at low 
T. We show that for the FFTXY and FFSQXY lattices ^nscha{T)~V{T) = 0{T^). 

We start with equ (0). Using the same notations as in Appendix A we have 

J{R) = I J^l cosa(^)e^^2/(-R) (Bl) 

X {6{R -e) + 6{R- e') + 6{R - r ) + 6{R -{f+e + e')) 
-6{R -e-f)- 6{R - e' - r)) 

where denotes the Kronecker delta symbol, where a{R) was defined in 

equs (|12|,|14|) and where the expression is written in such a way as to preserve the 
symmetry under the transformation e ^ e' . 

We now Fourier transform equ 



J<i =Y.\Jr\ cosa{R)e--^y^^^e-"^-^ (B2) 



R 



setting 



R f e,e' 

x{6{R - e) + 6{R - e') + 6{R - r) + 6{R - {f + e + e')) 
-6{R -e-r)-5{R-E' ~ r))e-''^-^ 



B = -e) + 6{R- e') +6{R-r)+ 6{R -{f+e + e')) 

R 

-5{R -£-r) -6{R-6' - f))e-"'-^ 
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we see that B = e-''^-^ + e''^-^' + e'^'^-^ ^ e''^-^^^^^^'^ - e-^^-^^+^ - e-^9-(^^+^^) 
Since we wish to compute 

r(T)= hm i(J(0) - J^.„-;) (B3) 

we may expand B to second order in g: 
5 = 2-(g.£)(g.£"^) + 0(g3) 

and thus : 

r(2^) = J E l^-H cosa(£-)(w;.£)2e-^^(^^ (B4) 



2 

'4T 



^g-i(y(£0+y(e^)+y(r+e>e^)+S/(r)-j/{r+e:)-j/{r+e')) 



Similarly, Jnscha{T) equ (|6D is given by 

iNSCHAiT) = W \M cosa{e){e.u-'fe--^y^'^ (B5) 



2 . 

E E(^-"'^^) (^•«'^) (I COS I I cos a(e') - J^-J^) 

We see that the difference between the expression of r(T) and the expression of 
iNSCHAiT) comes from the term proportionnal to cosa{e) cos «(£'). 
With the notations: 



2 

e 

r e.e' 
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^{y{e)+y{£')+y{r+£+S')+y{r)-y{r+s)-y{f+g')) 



xe 2 

y^Q-\{yi.^+yi.^)+yif+i+^)+y{^-y{r+e)-y{f+^)) 



1nscha{T) and T{T) read : 

1nscha{T) = A-^C + 
V{T) = A + 

Expanding A, C and S m.T yields: 

C = C° + C^T + C^T^ + and S ^ + S^T + S^T^ + 0{T^) 

We set g{r) = y{r)/T such that g{r) approaches a finite hmit as T — > 0. 

to order T° : 

C° = I ^ ^ I Je-I COS I I COS (ff.Ma;) .Ux) = 

r e,e' 

is clearly zero when one sums over e and ^ 
5° = for the SQ lattice because w = -1 so that Y.fW''-^''y = 0; 5° = for the TR 
lattice because Y.sJsi^-'^x) — 



to order T 



1 . 



.]^{g{e) + g{^) + g{r + e + ^) + g{r) - g{r + e) - g{r + ^))) 
= 
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{~){g{s) + g{^) + g{f+ e + ^) + ^(r) - ^(f + s) - g{f+ ^)) 
= 

Both and equal zero, owing to the parity in e and and using the fact 
that Y,fg{f+ a) = J2rg{'>^ whenever a is any vector connecting sites of the lattice. 

to order : 

= ^Xim cosa{e)\Jp\ cosa{e'){e.Ux){e* .u^) 

.^{g{e) + gie') + g{r + e + s') + g{f) - g{f+ e) - g{f+ e')^ 

.^{g{e) + g{^) + g{f+ e + ^) + g{f^ - g{r + e) - g{r + e'))^ 
Using the same properties as for the terms of order T we find: 

.]^g{r)g{f + £ + ^) 
Also expanding A to order T gives : 
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From these calculations we deduce that 'jNSCHAiT) = T(T) to order T { C = 
0{T^) ). At this order we may simply replace by in the expression of g{r) and 
we find that (equ (|28|) ) 

T 

1nscha{T) = 7o(1 - — ) (B6) 

with 

• For the triangular lattice 70 = v^/2J and T^o = — ^-j^fJ ~ 1.975J 

• For the square lattice 70 = V2/2J and T^o = ~ 2. 075 J 

~2 '2n 
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